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TITLE OF THE INVENTION 
DFT LEAKAGE REMOVAL FOR NON-COHERENTLY SAMPLED SIGNALS 

CROSS-REFERENCES TO RELATED APPLICATIONS: Not Applicable. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 
DEVELOPMENT: Not Applicable. 

Reference to Microfiche Appendix: Not Applicable 
BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

This invention relates generally to methods and equipment for testing electronic 
devices, and more particularly to techniques for accurately analyzing the spectra of non- 
coherently sampled test signals. 

DESCRIPTION OF RELATED ART INCLUDING INFORMATION 
DISCLOSED UNDER 37 C.F.R. 1.97 AND 1.98 

Test programs for automatic test equipment (ATE) commonly measure the spectra 
of signals sampled from devices under test (DUTs). Fig. 1 shows a conventional test 
scenario. An automatic test system 1 10 includes a host computer 1 12 that runs a test 
program. The test program activates a signal source 1 14 to apply a stimulus, which 
generally comprises one or more test tones, to the input of a DUT 120. A capture 
instrument 116 samples an output signal from the DUT 120 as it responds to the stimulus. 
The output signal from the DUT 120 includes the test tones, as well as noise and 
distortion. Tester software computes a Discrete Fourier Transform (DFT) of the output 
signal to produce a power spectrum. The magnitudes of the test tones are reported as the 
power levels of the DFT at frequencies corresponding to the test tones. Noise and 
distortion can be reported as the power levels of the DFT at all other frequencies. 

As is known, "leakage" error manifests itself in the DFT of a sampled test signal 
whenever the sample clock is not "coherent" with the tones of the input signal. Leakage 
is the mathematical consequence of performing a DFT on truncated tones — i.e., tones that 
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do not complete an integer number of cycles within the sample window. Leakage can be 
observed as an erroneous broadening of spectral lines, a creation of false peaks and 
troughs (lobes), and a general elevation of the power spectrum's noise floor. A sample 
clock is "coherent" with a test signal if its period multiplied by the number of captured 
samples is a precise integer multiple of the period of each tone in the test signal. 

An important characteristic of ATE is the ability to measure accurately noise and 
distortion from devices under test. Leakage directly impairs a test system's ability to 
measure noise and distortion by elevating the apparent noise floor of a DFT to a point 
where the noise and distortion from the device can no longer be observed. Leakage also 
impairs the test system' s ability to measure the magnitudes of the test tones themselves, 
because leakage induced from one test tone affects the magnitude of that test tone, as well 
as the magnitudes of all the other test tones. 

Several techniques have been used to reduce leakage. One technique is to 
multiply the sampled sequence by a windowing function that gradually tapers the 
sampled sequence to zero at its endpoints. The windowing function forces the windowed 
sequence to be periodic, and therefore coherent, within the sample window. It 
accomplishes this, however, at the expense of distorting the spectrum of the sampled 
signal and increasing the number of samples that must be captured, and therefore 
increasing the test time. 

Another technique is to mathematically convert the actual sampling rate to a rate 
that is coherent with the tones in the sampled signal. Sample rate conversion works by 
interpolating between actual points sampled at one rate to mathematically construct a 
series of points that appear to have been sampled at another rate. Although sample rate 
conversion can reduce leakage, it requires significant computation time and its accuracy 
can suffer from interpolation errors. 

Still another technique is to vary the rate of the actual sampling clock to ensure 
coherency. For example, in Fig. 1, values Ki and K 2 of dividers 124 and 126 are 
programmed to ensure that the width of the sampling window precisely equals an integer 
multiple of the period of every tone in the sampled signal. This technique is effective but 
requires expensive and complex hardware. It is particularly expensive when a tester 
includes a large number of sample clocks, as is often the case. 

Manufacturers of automatic test equipment seek to improve their products by 
providing less costly solutions to conventional testing problems. Great benefits can be 
gleaned from increasing tester performance while decreasing tester cost and complexity. 
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To this end, there is a strong incentive for inexpensively and conveniently reducing 
leakage to allow test systems to accurately measure the spectra of captured test signals. 
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BRIEF SUMMARY OF THE INVENTION 



With the foregoing background in mind, it is an object of the invention to measure 
accurately spectral components, such as noise and distortion, by reducing the effects of 
leakage on the spectra of non-coherently sampled test signals. 

To achieve the foregoing object, as well as other objectives and advantages, a 
technique according to the invention includes modeling the spectral components of one or 
more expected tones of a test signal, including the effects of leakage, based upon the 
frequency of the expected tone or tones a plurality of known sampling parameters. A 
DFT is taken of the sampled test signal, and the DFT is adjusted based on the modeled 
spectral components. The adjusted DFT is substantially leakage-free and directly reveals 
the components in the test signal, including low-power components that would ordinarily 
be lost in the leakage errors. 
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BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 



Additional objects, advantages, and novel features of the invention will become 
apparent from a consideration of the ensuing description and drawings, in which — 

Fig. 1 is a highly simplified block diagram of an automatic test system with which 
the invention may be employed; 

Figs. 2A-2D are spectra of different test signals shown to illustrate the principles 
of the invention; 

Fig. 3 is a simplified flow chart showing a process according to the invention for 
measuring the components in a non-coherently sampled signal; and 

Figs 4A and 4B are charts showing the performance of the process of Fig. 3 for 
measuring signal components substantially without the effects of leakage. 
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DETAILED DESCRIPTION OF THE INVENTION 



The technique according to the invention virtually eliminates leakage from the 
spectra of non-coherently sampled waveforms when the error in the actual sampling rate 
5 is known. To introduce the principles behind the invention, consider first a discrete-time 
test signal x[k] having a continuous-time spectrum X(f). Consider also a discrete-time, 
rectangular window function, w[k], which is defined as — 

w[k] =l,k = 0toN-l 
= 0, otherwise. 

10 As is known, the continuous-frequency spectrum of w[k] is given by — 

W(f) = e <--> 2 «^- m ) sin(oifN) / sin(7uf), (1) 



13 which is the discrete-time version of the well-known sin(x)/x, or sine, function. 

Analyzing x[k] over an interval of N samples effectively multiplies x[k] by w[k] in the 
discrete-time domain, which corresponds to a convolution of X(f) with W(f) in the 
15 frequency domain. 

For illustrative purposes, we assume that x[k] is a single pure tone A$w(27ikF t ), 
where the magnitude A is a complex number (thereby allowing x[k] to have a non-zero 
phase). The spectrum of x[k] has two frequency components: a first component at F t 
having amplitude A/2 and a second component at -F t having amplitude A*/2 ("*" denotes 
20 complex conjugate). The convolution of X(f) with W(f) therefore produces a pair of 
frequency-shifted rectangular window functions, which add together to yield the 
spectrum of the sampled tone. This spectrum can be expressed in the following general 
form, 

Xu{f) = X(f) ® W(f) = A/2 (W(f- F t )) + A*/2 (W(f - (F s - F,))), (2) 

25 where F s is the sampling rate and "®"designates convolution. By dividing each 
frequency by F s , we can normalize (2) as follows: 

Xwifv) = A/2 (W(fn - Ft/F s )) + A*/2 (^(f N - (1 - F t /F s ))). (3) 
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Fig. 2A shows the magnitude of Xw(f) as a function of normalized frequency, 
where F s = 1, and F t = 3/16. As shown in Fig. 2A, A^f N ) has a first peak at 3/16, which 
corresponds to F t , and a second peak at 13/16 (1 - 3/16), which corresponds to F s -F t . 

Fig. 2A might appear odd at first glance, because it shows the "continuous" 
spectrum 210 of the discrete-time signal x[k]. This is seldom observed. Instead, we 
normally observe what amounts to a sampled version of JG^fk) in the form of a DFT. As 
is known, a DFT is a "discretized" version of the continuous Fourier transform, which 
effectively "samples" the Fourier transform at regular intervals. The points 212 in Fig. 
2A represent the DFT of the sampled test signal. When using a DFT, we see only these 
points 212, and the underlying structure 210 of the Fourier transform is masked. It is 
crucial to recognize, however, that the underlying structure 210 is the actual spectrum of 
the windowed, sampled tone. The DFT is merely a group of regularly spaced samples 
taken from the larger structure. 

Xw(f) has the property that it peaks at precisely the tone frequency (and its 
reflection from F s ) and goes to zero at all other points where fN is an integer (See Eq. 1). 
When x[k] is coherently sampled, as it is in Fig. 2A, these peaks and zeros fall at 
precisely the points where a DFT would sample the spectrum. The resulting DFT returns 
non-zero values only at the tone frequencies, and returns zeroes at all other frequencies. 

This simple pattern breaks down, however, when x[k] is sampled with a non- 
coherent clock. In Fig. 2B, the normalized frequency of x[k] is changed from 3/16 to 
3.2/16, resulting in a non-coherent capture. Although the Fourier transform accurately 
reflects the magnitude of the sampled signal's spectrum, its peaks and zeroes do not 
coincide with the samples of the DFT. Rather than having the simple form that it did in 
Fig. 2A, the DFT of Fig. 2B suffers from leakage. 

Ordinarily, one does not know in advance the frequencies of signals upon which 
DFT's are taken, and leakage errors impair DFT accuracy. In ATE, however, the reverse 
is generally true. Test frequencies are almost always known in advance. A test engineer 
applies known stimuli to a DUT and samples the DUT's output at a known sampling rate. 
Output signals from the DUT generally have the same frequencies as input signals, or 
related frequencies such as harmonics and mixing products. Conventionally, test 
engineers have applied this knowledge of test parameters to set up coherent captures. But 
because ATE presents a case in which frequencies are known in advance, it also presents 
an opportunity to correct for leakage. 
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Fig. 3 shows an overview of a process according to the invention for correcting 
leakage errors in the DFT's of non-coherently sampled test signals. At step 3 10, a test 
program directs a test system to apply a stimulus to a DUT and sample a response from 
the DUT. The response includes at least one tone having known frequency. At step 3 12, 
the test program runs a DFT on the sampled response. At step 3 14, the test program 
models the spectral components of the known tone or tones, including the effects of 
leakage induced by the known tone or tones. The modeled components are then applied 
to corresponding components of the DFT, at step 3 16, to yield an adjusted DFT. The 
adjusted DFT provides a substantially leakage-free spectrum of the sampled test signal, 
which the test program can evaluate (at step 3 18) to determine whether the DUT passes 
or fails. 



Single Tone Correction 

Applying this process to a single-tone test signal, we first consider the 
requirements for a coherent capture. To sample a test signal coherently, sampling 
parameters M (# of cycles of input tone), N (# of samples), and F s (sampling frequency) 
are established such that M and N are mutually prime, and F t /M = Fg/N, or equivalently, 

F t = MFs/N. (4) 

Using these values, the DFT bins are spaced at integer multiples of Fg/N. Now assume 
that the capture is non-coherent, and that the sampling frequency is not exactly F s but 
rather F s ': 

F s ' = F s /(l+a), (5) 

where a is the ratio error. Rearranging terms and substituting F s from (5) into (4), we 
get— 

Ft = (l+ct)F s 'M/N. (6) 

We can express (6) in normalized form by dividing all frequencies by F s ': 

Ft/F s ' = (l+a)M/N, (6) 

where M/N is the normalized signal frequency. When F t /F s ' from (6) is substituted into 
(3), we obtain the following expression for the Fourier spectrum of the sampled, single- 
tone signal: 
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JGKfc) = A/2 (Wtfii - (l+a)M/N)) + A*/2 (Pf(f N - (1 - (1+aJMZN))). (7) 

Equation 7 can be rendered in a form analogous to a DFT by "discretizing" Xw($i), i.e., 
by sampling it at frequencies f N = 0, 1/N, 2/N, k/N, (N-l)/N, as follows: 

M = A/2 (^(k/N - (l+cc)M/N)) 

+ A*/2(^(k/N-(l-(l+a)M/N))) (8) 

Equation (8) thus represents the discrete spectrum, including leakage, of the expected 
sampled test signal. 

Although the frequency of the expected test tone is precisely known in advance, 
its amplitude is generally not known (a DUT generally does not have unity gain, so even 
if the amplitude of the stimulus is known, the amplitude of the response is generally not). 
To model the amplitude of Xw [k], we apply an actual value taken from a DFT obtained 
using substantially the same values of M, N, and F s ' that we used in modeling (8). We 
apply the measured value of the DFT in bin M (the tone frequency) to (8), as follows: 

X w [M] = DFT [M] = A/2 (fT(M/N - (1+O0M/N)) (9) 

+ A*/2 (W(Mm - (1- (l+ot)M/N))) 

We use bin M because it contains the largest expected tone, and thus can be measured 
more accurately than bins containing tones with much smaller magnitudes, which are 
more susceptible to noise and other errors. Other bins could theoretically be used, but 
with less satisfactory results. We can now solve for the two unknowns (A and A ) to 
complete our representation of X w [k]. It may appear at first that (9) prescribes one 
equation in two unknowns. But by separating (9) into real and imaginary parts, we obtain 
two equations in two unknowns: 

Re (DFT [M]) = Re (A/2 (W(M/N - (l+a)M/N)) (10) 

+ A*/2 (W(M/N - (1- (l+a)M/N))) ) 

Im(Z)Fr[M]) = Im(A/2(^(M^-(l+a)MM)) (11) 

+ A*/2 (W(WN - (1- (l+a)M/N))) ) 

We can simplify the notation as follows: 

Let E A = W(MM - (l+a)M/N)/2 
Let E B = W(MIN - (1- (l+a)M/N)))/2, 
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which allows us to rewrite (10) and (1 1): 

Re (DFT [M]) = Re (AE A + A*E B ) 
Im (DFT [M]) = Im (AE A + A*E B ). 

We can also rewrite the complex products as follows: 

AE A = [Re(A)Re(E A ) - Im(A)Im(E A )] + j [Im(A)Re(E A ) + Re(A)Im(E A )] 
AE B - [Re(A)Re(E B ) - Im(A)Im(E B )] + j [Im(A)Re(E B ) + Re(A)Im(E B )] 
A*E A = [Re(A*)Re(E A ) - Im(A*)Im(E A )] + j [Im(A*)Re(E A ) + Re(A*)Im(E A )] 
A*E B = [Re(A*)Re(E B ) - Im(A*)Im(E B )] + j [Im(A*)Re(E B ) + Re(A*)Im(E B )] 

By further recognizing that Re(A*) = Re(A) and Im(A*) = - Im(A), we can express (10) 
and (1 1) as two equations in two unknowns — Re(A) and Im(A): 

Re (DFT [M]) = Re(A) [Re(E A ) + Re(E B )] - Im(A) [Im(E A ) - Im(E B )] (12) 
Im (DFT [M]) = Re(A) [Im(E A ) + Im(E B )] + Im(A) [Re(E A ) - Re(E B )] (13) 

Because all the terms in (12) and (13) are real numbers, these equations can be solved 
directly using linear algebra. 

Once Re(A) and Im(A) are known, we can solve for A and A* by recognizing that 
A = Re(A) + j Im(A) and A* = Re(A) - j Im(A). These values of A and A* can be applied 
to (8) to generate any of the values of X w [k]. 

Until now we have assumed that the sampled signal is a pure tone. As we know, 
however, real test signals also include noise, distortion, and possibly other low-power 
components (herein referred to collectively as "NAD"). We can measure NAD by 
recognizing that Xy^\ should precisely match the DFT, except for the NAD in the DFT. 
Mathematically, we can express this as follows: 

NAD[k] - DFT[k] -X»{k], k * M, (14) 

where NAD[k] represents noise, distortion, and other low-power components in bin k, 
and DFT[k] represents the value of the DFT in bin k. Because Xy{k] includes the 
predicted leakage error, NAD[k] is substantially free of leakage and accurately represents 
the noise and distortion of the test signal. 
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Multitone Correction 

This approach can be extended to test signals comprising any number of tones. 
Assume that the test signal contains p different pure tones, wherein each i th of the p tones 
has an amplitude of Ai, a normalized frequency of Mj/N, and a ratio error a\ (ratio errors 

5 for different tones are generally the same, but may differ, particularly in microwave 
testing, and so are allowed to differ here). In the frequency domain, each i th tone 
corresponds to a pair of frequency-shifted rectangular window functions — a first at Mj/N 
and a second at 1 - Mi/N (see (8) above). At the bin corresponding to each i th of the p 
tones, the modeled discrete spectrum Xv{k] equals the sum of the contributions of 

10 rectangular window functions for each of the p tones: 

J^[k]= S [Ai/2(»tk^-(l+ai)Mi/N)) (15) 
+ Ai*/2 (W(k/N - (1- (1+aOMi/N)))]. 



Q 



Fig. 2D shows a graphical example of this summation of window functions for 
two tones (in particular, a summation of the tone in Fig. 2B with the tone in Fig. 2C). As 
p 15 with (8) above, (15) does not reflect the amplitudes of the test tones. As we did above, we 
apply data fropi a DFT of the sampled test signal, acquired using substantially the same 
values of M^N^ao^BUhat we used in modeling (15). We then apply the value of the 
DFT in each bin Mk (each tone frequency) to (15), as follows: 
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X w [M k ] = DFT [M k ] - I [^12 (W(M\JN - (l+aOH/N)) 

i=1 

20 (16) 

+ Ai*/2 (W(M k fN - (1- (l+ai)Mi/N)))]. 

Equation (16) is the multitone version of (9), and degenerates to (9) for p=l. 
Because there are p different tones M k , (16) expands into p different equations, one for 
each of Mi, M 2 , . . M p . Each instance of (16) can be expressed, as in (10) and (1 1), as 
25 two equations, one for Re (DFT [M k ]) and one for Im (DFT [M k ]). This expansion yields 
2p equations in 2p unknowns — Re(Ai) and Im(Ai). Each term in the expansion is a real 
number, which allows the system of equations to be solved using conventional linear 
algebra. 

As before, we can determine each complex amplitude Ai by recognizing that Ai = 
30 Re (Ai) + j Im (Ai) and Ai* = Re (Ai) - j Im (Ai). Once we know Ai and A;*, we can use 
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(15) to solve X w [k] for any value of k. We can then solve for each of the NAD 
components: 

NAD[k] = DFT[k] -X^k], (17) 

where 0 < k < N/2, and k c£ {Mi, M 2 , . . ., M p }. Equation 17 can be used for determining 
5 the leakage-corrected frequency component of any bin of the DFT. 

So far we have considered only leakage errors caused by expected tones in the 
sampled signal. One might question whether the low-power components themselves 
contribute leakage errors. Although all of the frequency components in a non-coherently 
sampled signal (even noise and distortion) contribute leakage, the amplitude of leakage 
10 errors from any component is significantly lower than the amplitude of the component 
itself. Leakage from low-power components such as noise and distortion is therefore 
well below the level of these components and can be safely ignored. 

In the preferred embodiment, the host computer 1 12 performs the mathematical 
operations under the direction of a test program. A DSP processor or special purpose 
15 hardware could be used in coordination with the host computer. Many of the 

mathematical computations required by this technique are identical for similar tests. 
Preferably, these identical computations are performed and stored for ready access by the 
^ tester during test execution. These terms may be stored when the test program is created, 

before the first execution of the test program, or the first time the test program is run. 

20 A simulation was performed to test the multiple-tone correction algorithm. A 5 12 

sample signal was created with 3 1 tones in bins 65, 66, . . . , 79, 8 1 , 82, . . . , 96. Each tone 
was generated at amplitude 1.0, and a random starting phase that was uniformly 
distributed from 0 to 2%. In place of the missing tone in bin 80, a spur was introduced at 
amplitude -100 dB to simulate intermodulation distortion. Uniformly distributed noise 
25 with peak-to-peak amplitude of 20 uV was added to each sample. Two sample sets were 
then generated: one with no coherence error, and the other with a coherence error (a) of 
1,000 ppm (parts per million). Fig. 4 A shows the smeared and corrected spectra of the 
capture. 

Three metrics were calculated to demonstrate the accuracy of the algorithm: 

30 1 . SNR: ratio of sum of signal bins to sum of all noise bins other than the spur in bin 

80; 

2. Spur: magnitude of spur in bin 80; 
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3 . SFDR: magnitude of largest noise spur, excluding bin 80. 

These metrics were calculated for the coherent data, smeared data, and corrected data. As 
can be seen in Fig. 4B, the coherent and corrected metrics are all within 0.1 dB of each 
other. Furthermore, the algorithm removed 66.9 dB, 60.5 dB, and 78.4 dB of coherence 
5 error from the SNR, spur, and SFDR metrics, respectively. 

Alternatives 

Having described one embodiment, numerous alternative embodiments or 
variations can be made. 

10 As described above, the technique requires that the tones in the input signal be 

known in advance. This requirement can be avoided, however, by scanning the DFT of 
M* the sampled test signal and examining the results. Peaks can be identified in the resulting 

2 power spectrum, and the frequencies corresponding to the peaks can be used as the tones 

fjj upon which the process is operated. The technique can also take into account attendant 

2 15 circumstances in identifying the tones of the test signal, such as harmonics of frequencies 

3 applied to the DUT and known characteristics of the DUT. 

As described above, the amplitudes of the tones of Xu{k] are established by 
M= applying one or more DFT measurements (See equations 9 and 16). Although DFT 

! H measurements provide accurate results, they are not strictly required. Alternatively, the 

U 20 amplitudes of the tones could be estimated rather than measured. Leakage errors would 
O- then be reduced in direct proportion to the accuracy of the estimated amplitudes. 

The technique has been described above as pertaining to the particular test 
scenario of Fig. 1. It should also be understood that the mathematical descriptions of the 
invention included herein can be expressed in various equivalent terms, via well-known 

25 mathematical transformations and/or substitutions, which accomplish the same or 

substantially similar results. Different expressions of the same underlying concepts are 
intended to fall within the scope of the invention. For example, we have defined a as a 
ratio error that represents the degree to which a capture is non-coherent. Non-coherency 
can be expressed in a number of different ways, however, the specific form of which is 

30 not critical to the invention. 

When a single tone test signal is involved, the technique for measuring NAD can 
be simplified as follows. By recognizing that leakage errors are exceedingly small near 
the peak of a tone's spectrum, we can substitute the peak actual DFT value for the single 
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tone case in place of the complex amplitudes A and A*. We can then determine the total 
noise and/or distortion as — 

NAD[k] * DFI\k] - [DFl\M\f2 (W(k/N - (l+a)M/N)) 
+ DF7IM]*/2 (W(kfN - (1- (l+a)M/N))) ], 

thus avoiding the need to solve for A and A*. This method of correcting leakage does not 
account for the leakage error of the tone bin of the DFT, however, and thus provides less 
accurate correction than the general solution described above. 

Each of these alternatives and variations, as well as others, has been contemplated 
by the inventors and is intended to fall within the scope of the instant invention. It should 
be understood, therefore, that the foregoing description is by way of example, and the 
invention should be limited only by the spirit and scope of the appended claims. 



14 



